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ABSTRACT 

The asymptotic solution for the transient analysis of a 
general nonlinear system in the neighborhood of the ' stability 
boundary is obtained by using the multiple-time-scaling 
asymptotic- expansion method. The nonlinearities are assumed 
to be of algebraic nature. Terms of order (where 6 

is the order of amplitude of the unknown) are included in the 
solution. The solution indicates that there always exists a 
limit -cycle. -i-The - limit: -cycle- is stable ^unstable) and exists- 
above . (below)— the stability boundary. if the nonlinear terms 
are stabilizing (destabilizing). Extension of the solution to 
include fifth order nonlinear terms is also presented. Com- 
parisons with harmonic balance and with multiple-time-scaling 
solution -of . panel f lutter-. equat ions -are- also.:! included? .■==. - 



ii 


FOREWARD 


This research was partially supported by a NASA Grant 
NGR-22-G04-030 . -Dr. -E. Car son_Yates , Jr., of NASA Langley 
Research Center acted as technical advisor. The author wishes 
to express his appreciation to Dr. Yates for the stimulating 
discussions made in connection with this work. The material 
covered in the main body of this report was presented as part 

of -a series— of- lectures given at -the - School.-.of Aerospace 

Eng ineer ingy - 7 of the Un i ver sity - of Rome ,”Xt a ly _*in Decembexyir. - - 


1972. 



TABLE OF CONTENTS 


PAGE NO. 


SECTION I: INTRODUCTION 1 

SECTION II: ONE -DEG REE -OF -FREEDOM SYSTEMS 7 

SECTION III: N-DEGREES-OF-FREEDOM SYSTEMS 13 

SECTION IV: FIFTH ORDER TERMS 23 

REFERENCES 34 

APPENDIX A: COMPARISON WITH HARMONIC BALANCE 

METHOD 3 6 

APPENDIX JB :7 “ SECOND-BORDER ^NONETNEARI-T-IES--^ 42 r: “ " 

APPENDIX^ C : THE .FUNCTION S-i (t-2) 46 



- 1 - 


SECTION I 
INTRODUCTION 

1.1 Regular Perturbation Methods 

This work deals with the stability analysis of nonlinear 
systems by a singular perturbation method, the multiple time 
scaling method. In order to motivate the use of this method, 
it is convenient to consider first the regular perturbation 
method . For the sake- of simplicity t h i s zme t hod —is _usedj..f or 
solving a simple equa tiohi^t-hetrBuf fln^equ abi-on^ 1 

X * X +£ X 5 =0 0<£^<1 ( 1 . 1 ) 

which- does- not -yield- any- instability, since it is conservative. 
The analysis of this equation* is useful- to under stand -the - 
meaning-of certain rterms " (secular termsi -andihence to mot-ivate--- 
the.-use-of -the— singular - perturbation' methods*^;.. As- well-known-,^- 
Eq.— (1.1) represent s-a - ma ss~ spring sy stem_(nni t mass with non- 

linear -spring) and - has - a -periodic - solution '(elliptic functions)-. 
It may be noted that Eq.~(l.l) has an energy -integral- - 

_X % ii 2 + £ X 1 = cons+ant (1.2) 

2 2 4 

which implies that the energy (and hence the: amplitude of vi- 
bration) remains bounded. 

According to a well known theorem ^-{Poincare— theorem--} , the. 
solution of Eq.' (1.1) depends analytical lylupon . the parameter 
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and hence may be expressed in the form 

x=x(t,e) = 2. <£" XnW 


X = x(.t,fc) 


- Z 

O 


g. Xr 


(-t) 


(1.3) 


where 


X n (-t) 


J_ O^K 

n L 


(1.4) 


The regular perturbation method consists in the use of Eq. (1.3) 
for solving Eq. (1.1). Combining Eqs. (1.1) and (1.3) yields 


f^X n +Z£"X n 

( X 0 + x 0 ) iH'VXtX, 3 ) <- £ i (^ 2 + x - +3X ° X '^-” 


= o 

(1.5) 


which implies 

X> x D = o 

(1.6a) 


x, +-X,.+ X c 3 = o - 

(1.6b). 


X z v X 2 X, = O 

(1.6c) 

The - solution of 

Eq. (1.6a) - — is: ... . 



X 0 = Q e Lt + 5e' ll: 

(1.7) 


where a is a complex member and a is its conjugate. Combining 


Eqs: (1.6b) and , Cl. 7) , one obtains 


X, 


+ x, + (a 3 e t31 * 3Q 2 ae lt + 3aa z e 


it , -srZ " it — 3 — l3t\ 

t a e J * O (i.8) 


The -solution of Eq . (1.8) .is 


X, = 


be 1 * +^Q 3 e l3t 1 1- Q-at e‘ 


t C ,T. 


(1.9) 


where C.T. represents the Conjugate Term of the term in brackets.- 
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Hence, the solution is given by 

X = x e +• £ x, fOC£ J ) = f(Q + 6 b) e Lt + £ £ Q 3 e l3t 

± £ l -£■ a* at e lt l + C.T, -»- 0 (.£. s ) (1.10 

The third term on the right hand side grows in time beyond any 
limit. Terms of this type are called secular terms . 1 Clearly, 
if the series is truncated to terms of second order , the so- 
lution does not have the expected properties since it is not 
periodic and does not remain- bounded in time. 


1.2 Singular Perturbation Methods 

In order to circumvent the above mentioned problem, several 
methods were introduced. These methods are called singular per- 
turbation methods and are based upon the concept of asymptotic 

2 3 

expansions. We 1 l^known rare-- the Krylov - and -Bol iubov. - method- - 

_ 4 _ -i. 

Caver ag-ing 'method-); -and- the , Cole, and xKevorkian . .method— { two- - - 

variable expansion) . More recent are the multiple scaling 
method (Kef s . 5-9) and"%he Lie - transform - method” (Ref s. “10-14) . 

The s e -las tr two ~~m e t hod s -ar e~u s ed- - i n R ef-.-J_ 5 to analyze- the stabi- * 
lity of nonlinear systems. The multiple time scaling is consi- 
dered here.* This method consists of assuming that the solution 
be of the form 


X - ) t £X, 

“ 2 - • • ) 


(l.ll) 


A comparison of the multiple time- scaling with the harmonic 
balance method is given in Appendix A. 
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where 

■+ = t n t (1.12) 

m 

are the multiple scales and are treated as independent variables. 

The functional dependence of the x upon the t is obtained by 

n m 

imposing that the solution contains no secular terms. It may 
be noted that Eg. (1.12) implies that 


dt ‘c Ot n d£ o dU 


(1.13) 


Combining Eqs. (1.1) , _ (1.11), and (1.13), one obtains 


t £ 


2 Jfdt, + • -•] (x 0 +£-X,j-efx,-) 


+ (,X 0 + £X, + £ Z X Z ) + £ X c + £ X, + . ... O 

which implies 


(1.14) 


- 3 ! x„ 

• + X = O - - 



(1. 

.15a) 

9 tV 






c) 2 x, 

t 2 3^ 

i- X, 

+ x* -o 

(1. 

.15b) 

3 -to 

3t 0 3t( 


U 

a 2 x. 

+ 7 9'x, + 

3 2 x 6 - 

:+ 2 & + x 2 +3*x r 

0(i. 

.15c) 

3 -to 

d4 o '0t ( 

s-tr 

c¥ a dt z 2 


9 4 


The solution of (1.15a) is 


^ lA 0 IT “ 

= Qe + a e 


where a and its conjugate a are functions of t.-,. t 
Combining Eqs. (1.15b) and (1.16) , one obtains 


(1.16) 


3 


— + ><(= - \?.i 2 a e dt 

di e 1 dt( 


3 *.3i 0 

+ ct e 


+ 3 crae it ]+c.T. d. 17 ) 
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where C.T. indicates the conjugate term. In order to avoid 
secular terms, the condition 

2 L §3 + 3cfa=o 

must be imposed. The solution of Eq. (1.18) is obtained by 

setting • 

Q = lale l(f 

and separating real and imaginary parts. This yields 


(1.18) 


or 


where 


D IqI, ■ _ q 
dU 


■ Q«e 


i- 2 


cl = IqI e 




(1.19) 

(1.20a) 

(1.20b) 

(1.21) 

( 1 . 22 ) 


is a function-of:-t 2 , _t 3 , If Eq. il. 18) is satisfied^ 


the solution of Eq. (1.17) is given by 

cK 1. ^ 1 


Vbe 1 ^ r gcre ~ + c.t.: 

Combining Eq. (1.11), (1.16), (1.21) and (1. 23 )-,- one obtains 

'•<.VfQ o a,,U 


(1.23) 


K = a 0 e 


£ (be H ° t- £ q» e^) + Die) + C.T. (1.24) 


which is equivalent to 

X= (Q 0 +£b + ...) e 1,at +• £ ( g° +•-..) e L5<jt -i- C.T. + 0(c J ) d/25) 
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with _ 

60 = 4 {1 ‘ 26) 

The analysis may be continued with the solution of Eq. (1.15). 

This yields another term in the expression for u) and a harmonic 
i5ot 

of the type e 


1.3 Comments 

It may be noted that the Taylor series of Eq. (1.26) yields 


X = Q 


,e Lt + 


b 


e lt +^Q 3 e l3t 


t 


if Qo 


Qote lt ] + 0(.t") 


(1.27) 


in agreement, with- Eq.“ (1.10). -Hence, - "the appearance-of- the - 
secular terms-does -not mean that the series solution does actually 
diverge if one considers the whole series , just in the same manner 
in which the whole series t - <t 3 /3!) + (t 5 /5!)... converges to 
s in- As- it- -is impossible -to bu i Id. -enough— terms_by successive 

approximations.- -to be able-to- ascer-ta-in=this : fact,-- the , appearance-^., 
of secular terms renders the method impracticable," (Ref. 1, 
p. 219)”. — The term te^"^ is due - to- the^f-act.xthat— the -difference . — 
between -two -sinusoidal-tune tions—having amplitude -and frequency - - 
slightly different, grows ( initially ) linear in t. Obviously, 
solutions which “contain secular terms are meaningful only if the 
complete series is known. Hence, they have limited interest 

from -the .practical point of view. - 

It may. be noted that the Duffin Equation has-been considered 

only as an example in order to introduce the concepts of secular 
terms - and singular perturbation methods.. — The r e suits -obtained . 
above- may be obtained with, more- elementary methods, -such- as 
Lindstedt's small perturbation method (Ref. 1, p. 224). 
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SECTION II 

ONE-DEGREE -OF -FREEDOM SYSTEMS 


2.1 Vander Pol Equation 

As it will be shown in Section III, a characteristic 
behavior of nonlinear systems in the neighborhood of the stabi- 
lity boundary is given by the existence of limit cycles. The 
simplest equation to yield a limit cycle solution is the Van 
der Pol equation. 

. 2 

Consider the -two -forms of the .Van— der-^Pol-equation 


CL K _ g: \ 

at 1 cu r c 3 


(&)' 


+■ X = o 


and 




(Z 

d£ 


( 2 . 1 ) 


( 2 . 2 ) 


Equa"tiori.i(2 . 2J : can- be- obtained . by differentiating Eq.v;(]2 . 1) 
with urespect ..to. t .and setting^Tr-j 

w - gU: 

' dbb 


(2.3) 


For ^simplicity , -the -multiple scaling- method is applied to Eq.._ - 

(2.1). Similar results are obtained for Eq. (2.2) (Ref. 5). 

Combining Eqs. (1.11), (1.13), and (2.1) yields 

— \ 2 . 2 . 

(dt* 2 M^dU + + + *-•-•) 


= 0 


which implies; 


9jS 0 

9 1. 


* **=o 


(2.4) 


(2.5a) 
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§4' + X, = '2 Jp°. +• Qxo - i (3xoV 

0-io difit, £>t 3 WJ 


p » « « 


The solution of Eq. (2.5a) is given by 

x « Qe lt ' +-Q e" ut ° 


where a and a are functions of t^, t 2 , 
(2.5b) and (2.6), one obtains 


(2.5b) 

( 2 . 6 ) 

Combining Eqs. 


!*< ^ =-2i|2e lto + iQe li (2 . 7 ) 


di 


The condition for avoiding secular terms is 

-2^9- + Q -CfQ = O- 


.*<P 


By setting - 

Q= Idle 

and separating real and imaginary parts, one obtains 


( 2 . 8 ) 


(2.9) 


original page ® 2 1 - /a I i I q [ 3 = o 

OF POOR QUAhiTx . 


or (p and' 


?) ( P. •• o 

(5 


. ^ <*« ' 

where k and are functions of t 2 , 

of Eq. (2.1) is given by 


(2.10a) 
(2.10 b): 

(2.12) 


Hence, the solution 




(2.13) 


The function 


^ \fuKer^. - 


us^e'- 


(2.14) 


Setting u= iVjl one. obtains 2t! =v.-2:j^ i3 


at, _ >a> 3^ 

or u= 1+ Ke**' which is equivalent to Eq. (2.12). 


3JaL = _ 


(|a|-|a| 3 )*-'u + 
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with 

S K = sgn(K) (2.15a) 

T = ii-jUl Ik} (2.15b) 

is plotted - 'in Fig. 1 which shows“cl early that the solution of — - 
Eq. (2.13) tends to a limit cycle 

X 0 (-t) « 2 cos (-UC&) ■ < 2 - 16) 

independent of the initial conditions (except for the phase 

angle 4o ) . 


2 . 2 — Gene rallied an— der— Pol- Equations 

A one-degree-of-f reedom system which has many characteristics 
of the N-degrees-of-f reedom systems (see Section'III) “is the 
general i zed ;;Van der Pol equation. 


x 4-x =■ -z& [(/SrX-^x) +4 r (lT e x - x 1 j 


(2.17) 


Comb in ing-- Eq s i "fl . H) y -(1.13) and _-(2 . 17.) " and.- separating-dif ferent 
powers of yields 

x — O- 

aji._ ^ 

O-ts (2 . 18a)_ 



t X, 



(2.18b) 


The'solution of Eq. (2.18) is given by Eq. - (2 . 6) f “ which'“when- 
combined with Eq. (2.18b) yields 


d\ 
2. 
: 0 


j- 4-X 

&tr 1 


= -2[t^ e lit - + ;U K (-taV% 


3 id a e‘ <0 ) ~ 3 t ( a 5 e" 4 - 3 da. e iio )J + C T - 


(2.19) 
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The condition for the elimination of the secular terms is 

J& + v ifii) Q- H- (tfe + (. 6i) Ol O. =0 2 

By using Eg. C2.9) and separating real and imaginary parts, one 


obtains 


+/% £ ial +-<T e lal 3 =o 

+(3 X +$ x \ a r=o 


The -"sodution of Eq. (2.21a) _ is 

lal = j /fe ~ 

\J 1 4- Kc? 

with k function of t.. Substituting into Eq. (2.21b) and in- 


tegrating yields 


<p4( 




with W 0 _■ function ^of t 2 1 V- finally, combining :Eqs; _(1 . 11) 

and---(1.13) with iEqs .- .42 . 63:,.. (2 . 22)- r ,-andu(2 . 23)-, one --obtains.^. 

X= lal 2 4^ Vn.\<x\ 4- cpj ( 2 . 2 - 

Qe 

with cO t = -fix 4- 5x • 

The behavior of the .solution depends upon the .signs-of. fio_ .andJSe. 

Sp - sgn (/3j 

s r= s 9 n Ue.) 

*See footnote to Eq. (2.12). 


(2.25a) 

(2.25b) 
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Note that Eq. (2.22) is equivalent to 


lal = r^'Sf, Sy 

with 

T s 2 /y^e I l&l 


(2.26) 


(2.27) 


The discussion of the solution behavior is divided into 
four cases, depending upon the signs of (?>£ and . 

Consider first the case analogous to the Van der Pol 
equation,- namely ^ 0 and-Y^i>_ 0 (destabilizing linear terms 

and stabilizing nonlinear ^termsH^- In this case,.' Eq.-^ (2 . 2 6) ~ 
yields 


lal 

I |5e/«el 


i - -P, J /[ +s K e 


-T 


(2.28) 


with* f^ (T). [plotted "in Fig^il. f In this ’case^thetsolution tendsr:-,- 
to a limit cycle given by 

X=2|qL (2.29) 


with 


I (X loO 51 l 
kO = \ +■ <£ 6O4 

JU. laL +-^6 


(2.30a) 

(2.30b) 

(2.30c) 


Next, consider the case with both— (linear and nonlinear) terms 
stabilizing ( ^ > 0 and ft's. > 0). In this case, Eq. (2.26) ' 
yields 


LaL =r tT ) = M 




US. e 


(2.31) 


The function- f 2 (T) (real only for < 0, 'T > 0) , is plotted in 
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Fig. 2, from which it is evident that the solution is uncondi- 
tionally stable. 

Third, consider the case with both (linear and nonlinear) 
terms destabilizing (. |3 R ^ 0 and HTg < 0) . In this case, Eg. 
(2.26) yields 


let I 

/Ski 2 







-T 


(2.32) 


The function f 3 (t) (real only for K< 0, T ^ 0) is shown in 
Fig. 3, from which -it is a pparent“tha tithe-: solution— is-unconddr=“ 
tionaliy - unstable ,- :and-.goes to inf inity—in - f inite time . — 

Finally consider the case with stabilizing linear terms 
C > 0) and destabilizing nonlinear terms ( < 0) . In 

this case, Eq. (2.26) yields 

The- -function f 4 (T) is-plotted— in-Fig. -4 , -f rom_which it is apparent - 

that there exists an unstable -limit -cycle-- given by --- 

X = zlaL.,, CLOS + 

with 

Cl) - \ 

4L = ^ fiP 0 

Equations._.C2 . 35.) are analogous to Eqs. (2.30) i The subscript - oo 
is used to indicate that the unstable limit cycle is reached as 


(2.35a) 

(2.35b) 

(2.35c) 


t 


— cxO • 
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SECTION III 

N-DEGREES-OF -FREEDOM SYSTEMS 


3.1 Introduction 

Consider an autonomous N-degrees-of-f reedom system of the 


type 


X. + APOx = -P(x) 


(3.1) 


where A is a matrix, function of a parameter A , and £ (x) is 
the vector of the nonlinear terms. Assume that the linear 
system 

x+A(Alx = o < 3 - 2 > 

is such that all the eigenvalues p^ of the matrix, A given by 

Deit(p 1+ AU)) =0 <3-3) 

are stable (negative real part) for A A c , while for A— A 0 
there exists“a 'paii: of eigenvalues p { -±tto and for • re 

— this pair of eigenvalues becomes unstable ( po s i t iv e_ r ea 1 par.t 
The value A 0 defines the stability limit. For A =A 0 / the 

solution of Eg. (3.2) is given by " '-"A 


X = uae + u,a e d.t. 


(3.4) 


where D.T. indicates the Damped Terms (eigenvalues with negative 
real part), u is the eigenvector relative to the eigenvalue lu) 
that is the nontrivial solution of the system, 

( iu) It A(/0 W =0 (3.5) 

The stability analysis of nonlinear systems considered here . ; 
deals with the study of Eq. (3.1) in the neighborhood of the 
value A 0 . 
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Assume that the matrix A is an analytic function of A 
in the neighborhood of A 0 . 


A = A„ *£Ux-K) + -A'fa-O*---- • • (3 - 6) - 

x 

and the nonlinear terms vector is an analytic function of x 

-P - {-fn] { bnp^j XpXtj + 2. Xp X, x r + 

Z d. Pire x P x r x s + x P x^XrXj x t +- • • • j (3 - 7) 

3.2 Formulation of Problem — — 

Assume for s imp licity_b that the., nonlinear terms ;are :only of 
odd order*, namely 


b 


h n 



= o 


In this case, if £ is the order of magnitude of X, the 
solution is of the type 


(3.8) 


X X i «L*-X#+-£ s X 7 v 


(3.9) 


where x t' X 3' X 5' 


are functions of ^0^2' ^4' 


(for. 


x ^ ^ 6 

the odd-order scales do not appear in the solution ) . Hence, 


dM " dt 4 ■■ d-L* - 

Furthermore, it is cbnvenient to set 

A = A + A z +■ £ 4 A^ + - * • • 

o *- 


where 

A,= ±l 


(3.10) 


C3.ll) 

(3.12) 


*If the even-order nonlinear terms are different from zero, 
the analysis is only formally more complicated, but the results 
are essentially the same (Appendix B) . 
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and A ^ to be determined (see Subsection 4.3). Combining Eq. 
(3.1) with Eqs. (3.6) through (3.11), one obtains 


1 


( A 0 + £*A 2 )(£X,t-t ! X s +£. ff x 5 + --** ) 


— £ 3 -f 3 + £ S f s -I- - - ■ . • 


with 


" \A t A! 


(3.13) 

(3.14a) :■ 

(3.14b) “ 


and 


^ Xi >p 

-F 3 2 Cnp^r ( X«,p Xi^^X s,rA+^ Xl^'Xs^X ) 


(3.15a) 


+ p^' G^-(X'.P X ^ X ' 
Equal ion— (-3 .13) - i-mpl ie s - - ■ 


,r Xi,v X s t)^ 


( 3 . 15b)__. 


— ' f A X - O 
a-t c o1 

is **«.“*. 


- A,x 

Ot^ 


2*1 


- Ay - AX 


(3.16a) 

(3.16b) 

(3.16c) 


3. 3 A 1 Third Order -Solution • 


The solution of Eqs. (3.16) including third order. terras. 
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is considered here. The effect of fifth order terms is consi- 
dered in Section 4. The solution to Eq. (3.16) is given by 


X, -U ae L ° -fU 


(3.17) 


In Eq. (3.17), the damped terms have been disregarded, since, 
as shown for instance in Ref . 6 , these terms do not have any 
effect on the solution as ~t — . 

By combining Eqs. (3.16b) and (3.17), one obtains 


ir+Ax, =£ <3) a 


I3u)t 


di 


with 


'O '3 


(.0 _ tuif 

q Cue. — 


t ----- 


U • ' ■ A 


3 

t 


(3.18 ) 


01 


U Q : e 


C-oT. 


r=izo., r i 


- tivp U ^ lA y 


■^ t0 = }Z G ^ U P u % u r+ W p u U r+ u F (A u r J 


(3.19a) 

(3.19b) 


P<p“_ ... ‘ D r -b • r l 

The- condition! f or_ the . elimination -of the! ^secular' terms : is ;that r 


the— component of -the— vector 





(3.20a) 


in the direction of the vector u. is equal to zero. 6 In other 

words, the inner product between the vector 2 3 ^ and the eigen- 

-T 

vector u a of the adjoint operator (A is the conjugate, of _ the 

n 

transpose of A) , defined by 

(-LU)I +-A T )ci^ = O (3.20b) 

g 

must be equal to zero or 



= o 


(3.21) 



-17- 


Assuming for simplicity that the vector u A is normalized by the 
condition 

(u Aj u) = uju - 1 (3 - 22) 

Equation (.3.21) is equivalent to the condition 


2&- 4 Acl +5cCa=o 

€>■ U ' 

|S-(u«,A 2 u)= u*A 2 u 


with 


(3.23a) 

(3.23b) 

(3.23c) 


* 3 

Equation '(3. 23a) "is similar to Eq.- (2 . 20) rand_ the solution, is 
given by 


a=lale 5,1 


(3.24) 


with 


lal=: 




4 4- X-e 




- (3 X - - 


(3.2 5a.): 


(3.25b) 


wher e=4^ ' and - k- are^nf uncteions-of ^ , : _t - ^ . . . and-:^ R , r A j* :: 
and y are , respectively, the real and imaginary parts of /S 
and ft . Finally, the solution of Eq. (3.18), under the condition 
expressed by Eq. (3.23) is given by 


1 3 cot. 


x 3 = p 3 (3) e l3u ^° + (b 


(t) \ t ulio 

U^p 3 j ) e 


(3.26) 


(3) 


where b is a function of t^r 1 4 .... (see. Section 4) and P 3 w; aiid 
p are the solutions of the systems. 
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['i3a3l+Ajp 3 lS) =-f 3 l V ( 

[iu)I+A c ] p 3 l ° = a*a -u|-A 2 ua < 

respectively. The vector P 3 ^ contains an arbitrary constant 
which can be eliminated by setting, for convenience, 

(u A) p 3 o1 ) = u* T p 3 c,) = o ( 


3 . 4 ~ Comments 

Note. that, the solution of Eq.‘ (3.16b) is important not as 
much for obtaining the solution for X 3 (tg). (see. Eq._3. 2 6 ),.— but . 
for “obtaining the "functional- dependence of- x-r .upon -t^ . ~ Eq. (3. 
contains the function b(t 2 ....) which can be obtained only by 
studying -Eq.r (3.16c) for - the .unknown x^ (t 0 ) Hence,- considers 
the solution -.which .-can -_be written.as. 

X -£ (uQ.e. lU>t + Ud e" 1 ?* ) + . ^ ( £3 ) < 3 

where:a is given-by Eq.-~ (3.24)*.— In Eq. -(3.29)“, the damped .-vv 

terms-=have been- disregarded . Furthermore., still, neglecting 

higher order terms, £ is related to A by 

(3 

with ' X z - ± { , or 

X,_= sgn (A.-XJ 

It should be noted that by definition of X , the linear terms 
are stabilizing (which implies - ^>0) for A<X 0 (X 2 = : -() and 


.27a) 

.27b) 

.28) 


26) 

29) 

30) 

31 ) 

32 ) 
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vice versa for 
one obtains 


A >A 0 . Combining Eqs. 


p = ^(u^A'u) 


(3.23a) and (3.14a) 


(3.33) 


Hence /3 R > 0 for -1, implies that the real part of 

(u , A 1 u) is negative. Hence, the discussion based on the re- 
suits of Section 2, reduces to two cases only, > 0 and 

< 0, respectively. In the first case, X R > 0 (stabilizing 

nonlinear terms) , the solution is unconditionally stable for 
A A 0 . (Fig. 2) while for A >X 0 there exists a stable limit 


c y c 1 e_" = (F_ig.;- _1 ) given by Ir.,- 



e 


C3.341 


with CfU given by Eq. (2.30c). This behavior is summarized 
in F_ig . 5 : where - the amplitude of the limit cycle 



(3.35) 


is plotted : and “the differ ent- trends- -of the -solution - (for -A > \ Q „ 
and A < A 0 ) are indicated. 

In the second case, 0 (destabilizing nonlinear terms), 

the solution is unconditionally- stabl e_~ (Fig . :3) for_ A >A 0 • 

while, for A A 0 , there exists an unstable -limit cycle (Fig. 4) 
given by 


= J(X-X D ) h iue +ne J 


(3.36) 


with given by Eq. (2.35c). This behavior is summarized in 
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Fig. 6 where the amplitude of the limit cycle 


h 


L 


c. 



(3.37) 


is plotted and the different trends of the solution (for A>A 0 
and X^A 0 ) are indicated. 


3.5 Application to Aeroelasticity 

The system of equations for panel flutter is of the type 



y + Gy i Ky ^X Ey. = o : 

(3.38) 

with 





(3.39a) 


<3 * { 9 ^] = 1 '^ , r -A r 'PH r yp y<i Y r 1 

(3.39b) 

Equation 

(3 . 38 ) is equivalent ."to the -system - 



y = z 

(3.40a) 


Z + G z + K y + A E y = q 

(3.40b) 

which is 

.of the = type-of Eq.. -'(3.1) 



X-ffl 

(3.41a) 


Htl 

(3.41b) 


and .finally 



O ! I 
K+XE* d+ixG 

L i ' 


A o .*A'(-A-X>.. . 


(3.42) 



with 


L^o 


I 

g‘ 


1 « 
< 


A 


0 « o 

1 eTg'i 


where 


g 0 =g*|x 6 g, gStt. s and K * = K +XoE ° 

Furthermore, the eigenvector u is given by 


u = <- -- , 

Iw J 

where v and w can be obtained from 


(3.43a) 


(3.43b) 


(3.44) 


iu> v--^-w -gc^ -r, 

tco W/ + (G 0 w + K 0 v) = o 

This implies that v is the eigenvector'- of the original "system 


(3.45a): 

(3.45b) 


[-u£l..-Ki -td G 0 i- K 0 ]v ~ O 


while W is- given -by Eq. — (3 . 45a)- : -and hence - 


(3.46) 


u 


V 

TaJV 


(3.47) 


Similarly., the eigenvector u A is given by 

u A T [/ tol + = o 

Setting. 



one obtains 

<cow T + v t Ko = q 


(3.48) 


(3.49) 


(3.50a) 
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. . . A.T 

or, eliminating w 

v T [~^I + *-loG 0 + K 0 ] = o 

In other words, the vector v T is the left-eigenvector of the 

*A 


(3.51) 


— T 

original systems and the vector uu is given by 


u T = l_v T (iiAl +Gj! v T J= b^>v T K 0 |v 1 J (3.52) 


- T 


Finally, by using the first expression for u A r Eq. 
e q u iva ien t- to . 

u J U Lv 1 (tujI+G;)' ’! v T J 

= V T (2iu)1+G)v = \ 

while Eqs . (3.23) are equivalent to - 

A'u = X t lw Tl lV T J [g irgr] j- 
= v T i Er- -Et ^ G ) v * ^ 


(3.22) is 


V 

diov 


(3.53) 


(3.54a) 


with 


(3.54b) 


K ={Zl r C np<lr (u p u t u r > u p u % u r v U^U,) 

in agreement with the results of Refs. 6 and 7. 


(3.55) 
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SECTION IV 
FIFTH ORDER TERMS 


4.1 Introduction 

In this ■ section., - the ef feet— of -the ; fifth order—noniinear-- 
terms is analyzed. Combining Eqs. (3.16c) , (3.17) , and (3.25) 



(4.1) 


Furthermore,- according to Eq . (3.15b), 



(4.3) 
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By combining Eqs. (4.2), (3.17), and (3.26), one obtains 

£5 = {is,"] ‘ { Z> c„ PV - [a'e 1 ^ UpU^ + aa(u P u^ + UpU^) 

1[ .is) iiuiio 1 . 1 .. ' w\ iu>t„ ] 

JL p 3r e + 1 b u r + p 3 ,r je -+ c.T.J 


+- CL Up 0 


-l 2 u?to- 


+ [(^Upi^UrUsUje 


c5a>t e 


+ 


a^ lOpU^UrUsU* +'UpiJ <j) au5U f + u^u r u s u t + . . . 

+ UpU^UrU^uJe 13 + a’a* (u P u^UfU^Ut *■ u p u 4 Li r u s u t 
+ . .:.,kUp u*u..u s u t ) : e iuy *° + c. t. ] 1 ■; - 


f (5) tStDto J* 

S 0 + 1 


0 f* ^ t3tot 0 rt LU^to _ __ 

e + -K e +- CT. 


(4.4) 


where 


5 | j?np<^r to - u p u^ s (. U> la p "t ,.p 3 .b aa(upU^ tUpU^)x . 

P5- ? 

lb Ur + pj'r^) + 5/" OpU| p 3 “jb 2 -Srvp^rst OrdMUpU^-UrUi-U* 


Up U c^.U r U«, li-t- 4" *•• - »*+■- Up U^_U 


?y^. 

..u*uj.j ... 


(4.5) 


Similar expressions- hold for_£^ 1?^ and_£ c — ^ 
noted that Eq. (3.27a) may be rewritten as 
^ < 5 ) 3 /\ C3) 

Ps = Q P 3 


It may be 


(4.6) 
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with P 3 ^ 3 ^ independent of a, while combining Eqs. (3.27b) and 
(3.23a) f one may write 

with p - ^ and p 3 independent of a. In summary, it is possible 

JL J 

to rewrite Eq. (4.5) as 



(4.8) 


where-f is given by Eq. ..(3 . 19b). .or , -according to Eq. (4.3) 


= |Z,'C,p r ( UpU^Ur *• UpU^U, t UpU^U 


■»] 


{2 c„ 

( p^r - 


pc^r Up U<^ U 




(4.9) 

On the \o th er; f hand the - explicit- expressions • for f 3 and •< 
(independent of a) can -be obtained— by combinihg--E^s i -- {4.5)-,- 
(4. 6) ,-and~(:4 . 7) . Finally, combining Eqs. (4.1) and (4.9) and 
noting. -that-rac cor ding -to Eqsv--(-3 . 23) and— (4.7) , it is possible 


to write 



d± z 


A 


_ (0 

P 3 Q 


+ 


P?dcx 




p"a 3 a 


(4.10) 


(with p 3 


( 1 ) 




independent -of a) ,- -one. obtains — • 


2*5 +- ^ X 5 - 




Scot, 


+z 


(5) 


e +■ 


u) ..UtXiL 


Z'e 


-t-C.T. 

(4.11) 
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where 

Z'g 1 =-[u ^ + A ; ub --f s cl ^(o‘b 4-2 adb) + *-* 

+ (A<ut Pj 10 + A 2 p J “ , -$ s l,, )a 
+ ( j^+A.f,."' “W + rfa'] 


(3)-, 

Similar expressions hold for Z^ and; Z^ . 
The 'non- secular- terms condition is 



or , according - to Eqs^-_ (3 . 22) , - (3 . 23a) , - and — (3 .-23b):— 




where -: - 


~ 5a +Sa , a -i-S'h 3 a' 


da 

D4: 4, 



The~: solution of Eq. (4.14). is discus sed-.in Subsec.tion _ _.4..-2 
Finally, if Eq. (4.14). is satisfied, _ the solution of Eq. 
is giyen by 


(4.12) 

(4.13) 

(4.14) 

(4.15a) 

(4.15b) 

(4.15c) 

. 11 ) 
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= Q C5) 
5 K5 


‘i5a>t 


p 5 ^'e' ^pre + 


(3) 

5 


t 3totc 


( pg 0 tc«)e l 


cot* 


-b C.T. 


(4.16) 


00. 


where p 5 is the solution of the algebraic system 


[iK^l+Ajpf =- 


(0. 


(4.17) 


4.2 The Function b (t^) 

Consider Eq. (4.14) using Eq. (2.9) 

«P_ , 


setting 


q = Iq.I o 1 ' 


b = b'e L ^ 

up 


and dividing by e Y , yields 


+^(lar b '^2ta |1 b') 

= "&ldl .*-S'lar + <5" Cal Id d§ ■ 

OU- Qi v 


(4 . 18 ) . -- 
(4.19) 


(4.20) -_~ 


The" imaginary: part _of- Eq i™(3* 23) i yields :4see* also._.Eq. 2 . 21b)_ . 

--fr-W (4 - 21) 

Combining this equation with Eqv — (4.20) yields™ 

^^>^| Q rbuiar(.bbb') 

= Slal ^'lQl%&"lal s -gg-ilalM. 

D4-a ( a o'>\ 
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Next, by separating real and imaginary parts, one obtains 




S.lol 4 -S'IqI t5; lal’-Stol 

d±4 


§j^ ^b^UoTbi +2^|alX 

= & x I a I la.l 3 + <5 r " io. I fe - \a\d$ , 

- 

Finally, by setting — 

b' = B p 5 ]qI 
R a-t i 

bt = B x lQl 

and using the real part of Eq.~(3.23) (see also Eq. 2.21a) 
and- its derivative ..with -respect— to., — 


^laPdiQl _ 


ai* 


one obtains 


dB, 


a-t* 


©Bx 


- SrIqI 1 pi* » S R IqI 5 - QI q I /dig 

b lQl/ak 


^ =&r ^xhr+6;iai 4 - 

It may be noted .that , according to Eq. . (3.25a) , 
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with 



\ ^ tKe 11 


" 4 +K„e 


/Se 


Ko=-<Wl 5 e 

T = 2 

which yields 

gal = -/2> e IqI 3 K 0 e T 


and 

QIq I = __l p _ 1 3 Q^oj p-x. . . •• 

a-t 4 - - 2 |UI o-t*- .. ' 


(4.28) 


(4.29) 

(4.30) 


(4.31) 


(4.32) 


Combining Eqs.. (4.27a), (4.28), (4.31)-, and (4.32) yields 


d - 

c>? 


r 



t dK° "r- 

a-u '• 



(4.33) 


and hence by. integrating,— one-obtains^-.- 


B e = 

(3c 




2 5£ 


X'/u 1 0.4 ■ 



2K 0 _ cMr 4 ] L 


Note that 

\ 

Jc + K^e 1 


dT = 


c 


e -1 - x jMc-t-K 0 e' t ) 
c c 


(4.34) 
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where is a function of 


The condition for non-secular terms* is 




(4.35) 


Combining Eqs. C4.34), (4.35), (4.30), (4.25), and (4.24a) yields 


b' = 


- 

2 &?-l 


^ J 


+ 1 Ip' ' •>$«• ^- la l)(S E | a! t- Y e | ct! ?)-i 


= B * o) lal +■ B r l1) 1 qI 3 + (b^Iq! -+ B^| Q | 3 )i^fa 


where as mentioned above, is a function of * 


(4.36) 
. , while 


R <0.. 
R 




+ 




(Sr. 


D ^ _ - - (3fc CM 
Ti Op 




(4.37a) 

(4.37b) 

(4.37c) 


R ^ = L_ C " 

bR ^ 

Similarly., by imposing -the condibion- (Appendix C, Eqs; - C. 13. and 
C.14) 


*The secular terms must -be considered -for "t — in- the case-of - 
stable limit cycle , and for“t— ^-°oin the case of unstable d-irait... 
cycle. 
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DCk 

34 


-- s, -s; & * s; ( - 1: l s.# * ^ ( fj) - * 


11 V Sfc 


fc- 


(4.38) 


one obtains (Appendix C, Eq. C.18) 


bt= B4 0 ) Iq| 4- 


< 4 - 39 > 


where B ^ is a function of t~ while (Appendix C, Eq. C.17) 


A 


Ki R W 

R 
R 


u u) 
bj ,=- - 




b,' 1 - > fcfe l! > 

Or 

Finally, combining. Eqs._ ._(4 .19) , (4.36)., and— (4.39) yields 

b = B to b +- -t- ( B u <x t B^Va) 2^ lal 


(4.40a) 


(4.40b) 


(4.40c) 


with 


^ o<) . , o CK) 

B ~ B|5 +■ l Pj 


(4.41) 


(4.42) 


4.3- The Function a (t^Y t^ . . . . ) 

The solution of Eq. (4.35) is 
C-t^. 


K c = K,e 


(4.43) 


with 


G 




(4.44) 


while the -solution of Eq. (4.35) is 

2?o = ^4 +( Pi 


(4.45) 



(4.46) 


Combining Eqs. 


(4.43), (4.45), 

a= lale ltp 


(3.24), and (4.28) yields 


with 


Iq'= 

(P = . cO ^ ,+ Sx |g[| + 
Kr. 


(4.47a) 


(4.47b) 


Finally ^, it is possible— to determine— the -arbitrary -constant _ 
A^ introduced— in Eq. “ (-3 . 11) -- -Note^-that— A 4 - appears-- only -in 
which'is given by Eq. (3.14b) r 

A,.-X 1 A'<-A'-X 1 A J tA" ■ l1 - 181 

and that appears only in 5 - , given-by Eq.— (4 .15a) ..,™. - 


& ■ XU;!?® A_M .A u .jp ". A, p’t 


CO 


(4.49) 


where--®5c) -is the -value -for 6 at ~A 4 - _0 and ~ 0 >l is given-by Eq. 
(3.23a). Hence, Eq. (4.44) may be rewritten as 


= ^ ) (Sr + 


(4.50) 


and whence the exponent in Eq.- (4.47a) is 



If <5 0 has the opposite sign as ft ( which is positive for ^ 0 

and negative for > 0) then the limit cycle switches from un- 
stable to stable (or vice versa) for 



2(3 a Xt + (T 0 

-A*. 


(4 



This value of 6 must correspond to the knee of 
as a function of X , and hence must satisfy the 


the curve 
condition 


OX 

dt\ - 





— O 


Combining Eqs^r:(:4 . 52 ) and- { 4 . 53 ) yields — 


(4 


^ Q ~ 0 x> (4 

2 /S/e 

Applications of these- results to the panel-flutter problem are 
given -:in Ref 8 , where the equation - 

x +[x+ ) u(x 2 + x 1 ) +9 (Xf-xf] X +-.X -O (4 

(with ^X <-- 0, }x > - 0, <r : 0) is also.^considered - This.equa- 

tion yields "two limit cycles— (one stable and -one unstable) 
given by (exact solution! )'~ 

X - X K COS (t + <&) ( k - '/ 2- ) (4 

where - X x are the'two' positive -roots of 


52 ) 


53) 


54) 


55) 


. 56 ) 
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APPENDIX A 

COMPARISON WITH HARMONIC BALANCE METHOD 
A.l Introduction 

In this Appendix, the multiple time scaling method is 
compared to the harmonic balance method. For simplicity, the 
analysis is limited to nonlinear terms of third order (see Eqs. 
3,1, 3.7, and 3.8). 

First,— the system 
N 

X n -+-Z Xp — ?npejrL2^ p "X cj^X r (A. 1) 


where—f— is the order'of magnitude -of -X n is analyzed -by the--v_ 
harmonic balance method. The resulting system of algebraic 
equat ions -i s = then solved "by an algebra icr.pertur bat ion methods.; .. 
Finally^ lit is shown that --the -.re suits ^.coincide ^with. the ones — 
obtained -by using -the-multiple ."time scaling method . 

Consider tEq. - (A.l)-. a According -to the .Harmonic . Balancer - r 
Method-r-lwhich is equivalent- to the - Galerkin- method*) 1 , the solution, 
of Eq.-<A.l) can be obtained by setting 


x„ = Zix^e^Vxfe'^) 


(A. 2) 


■^fpi 


(where X^ is a complex number and • X n r is its complex 
conjugate) and- by “balancing” the firstrP harmonics of -the re-,- 
suiting expression:, i.e., by settingzthe coefficients of “ 
equal to zero for p = 1,2 . . . .P. For. P = 1, and. setting X n -A n 
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one obtains 

^ tc) +■ ^npXp “ ^npi^r fXpXcX r t X p X<^X r 

^ (A. 3) 

4- XpXcjXj + 0(t s ) 

and the complex conjugate equation. This system of algebraic 
equations can be solved by a standard numerical technique such 
as the Newton-Raphson method, for instance, (Refs. 7 and 14). 
Here, Eq. (A. 3) is solved by an algebraic perturbation method. 


A. 2 Algebraic .Perturbation Method . . 

Since y is of order S~ , set 
n 

X n = t x» t , + £ 3 Xn, 5 + 0U £ ) (A. 4) - 

By combining Eqs.,- (A. 3), -and— (-A. 4) and neglect ing_higher_order:_-_ 
te rm s-, -o ne- obtain s- 


a/ 


(A. 5)— 


The "’solution of Eq.~(A.5) is discussed in Section III _ (see :-Eq.r_-. 
3.5) which - shows that' a- solution- exists for _X= Xv and_.<u> a) 0 - 
and is given by 


where is the eigenvector of the matrix 
to the eigenvalue . Next set 


lA„ P (x Si 


£<)=■. uy 0 + e z ui z + 0(A 4 ) ' 
A - A 0 t A z t 0(e 4 ) • 


(A. 6) 
relative 


(A. 7) 
(A. 8) 


with X 2 =±4 (see Eqs. 3.11 and 3.12). 
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By 


using Eg. (A. 8), one obtains {see Eqs. 3.6 and 3.14a), 


Anp^ = AnplXt.) + fefr ) £* A* + 

X* Xn 


«/C 


u) 

p 


(A. 9) 


wi th 


Anp " A* p (Xj 


(A.-10)- 


Combining Eqs. _(A.3), (A. 4), (A. 7-)-, -and (A. 9) and— separating 

terms of. the same order of magnitude yields. 


Xn,< i- Xl j A n p Xp,{ O 


t ^o Xn;V + -2 A n p 


(A. 11a) 


(A.ll-b) 


with 


Z h = ”4 Xru ~ ^ A h p } Xpj 


P^ 1 


v Xp ( < Xr,< t X<^,1 X.r.1 +■ .Xp.(;X^ ( , X, 


J 


(A. 12) 


The solution. o.f Eq. (A. 11a) is given by Eq. (A. 6). Combining 
Eqs.- (A. 6) and (A. 12) yields- 


Z n = - tcd-i CL U n ~ <2.2 Anp U 


+• & & 2 ^Vp<^r ( Up Uc^ Ur Up U(^ Ur Up Ur ) 


(A. 13) 



Equation (A. lib) has a determinant equal to zero. Hence, the 
solution exists only if t. is orthogonal to the left eigen- 


vector 




(see Eq. 3.20b) of the matrix 


[A„n 


2 Ur i-t\ " “ ^ ^*0-2 j Vn Q- 2j ^ Anp 

+ CL (X 2j ^ ^ ^ ^ U(^ Uh) O 


or 


L CjD z XL - + /S> QJr -h S Ou CL = - Q 


(A. 14) 


(A. 15) 


with 


/3- 2r^A, P Up 

— 7 (A . 16 ) _ 

Z'^Un =, 

^ ^.p^rl.Up A ^ ^4 ^ \ ^ A 17 y . 

^-/»c — 

in agreement with Eqs.-43.22), (3 . 23a^-a-nd -(3 . 23 b) . Separating 

the real and imaginary parts, Eq. (A. 15) yields 

fitCL f ^T a a J =0 (A. 18) 

co z a. a. +■ tf z a 3= o (A.i9) 


Equation (A. 18) yields 
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( A* , Eq. (3.12) must be chosen so that a is real) while Eq. 

(A. 19) yields 

td z = -/3 X - ^ d 1 = -/5 X + & A (A. 21) 

tte. 

in agreement with Eqs . (2.30a) and (3.25b) , respectively. 

In conclusion, the harmonic balance method (combined with 
an algebraic perturbation technique) yields the same results— as 
the multiple time scaling. It may be noted that the results 
are more easilyiderived with, the 4 harmonic bal.an.ce.__ method than 

with' the : multiple --timet sea iingtme thodTir_-_However , :i the” harmonic-- 

balance- method is limited -to steady . state= 41 imit-cycle) solutions 
while the multiple time scaling method yields (with no additional 
computational. efforti--- also,: the envelope' of the transient res- 
ponse and, in particular, the stability (or instability) or the 
limitr cycle. - 

Finally 7 it may be worth noting a computational advantage 
tha t“ the-a 1 g ebr ale r per turba t ion . harmonic - ba la nc e . me thod _ { a s 
we Ilia s the multiple time scaling method)- has with respect to 
the Newton-Raphson harmonic balance method. Solving the linear 
system, - ~Eq. - (A. 5) , requires -less :time than -solving the nonlinear 
system, - Eq. (A. 3) . . Moreover.., . the nonlinear- system must be 
solved for different values of A" , while the linear system must 
be solved only once: the additional time to obtain the coeffi- 

cients and_.fr... is negligible. Hence , -the total time required 
by the’Newton-Raphson- harmonic balance-method is N- times, larger 
than -the one required by the algebraic-perturbation harmonic 
balance method (and multiple time scaling method) where N is 
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the number of points used to describe the curve which gives the 
amplitude of the limit cycle versus A . Note that the higher 
accuracy obtained by solving exactly Eq. (A. 3) (Newton-Raphson 
harmonic balance) is only apparent since the difference between 
the two results is of the same order of magnitude,- c , as the 
error in the original equation, Eq. (A.l). 



APPENDIX B 


SECOND ORDER NONLINEARITIES 


B.l Introduction 

In the* main- body of this- 'repor-t-JEq 1) is analyzed 
under the assumption that the even order terms are equal. to 
zero, Eq. (3.8). The solution of Eq. (3.16) including third 
order terms is considered in this Appendix without the assump 


tions b„ p ^= O and = O .The solution is 

assumed to be given- by 

X =:-£ X, (£«) ■ 

where--x^- / x 2 /'^ - • . ..are.-: functions r_©f---6g ,t 1> /t 2 .- Hence— - 


(B. 


d _ £L + a 0- + & '<±- 

dH ' dt 0 + £ dii S &K 


’Z' ■ ■ o • 


(B. 


_Q>mbinirig=.;Eqs Lr..{3 . 11 ,1(3. 6.1, 2(3 . 7.)- , - {B . I)-, -and B. 2 )_, -One 
obtains "z- _ 

( kK&t c k ■ 

+ ( Ad^.?A 2 M.y- ■ )(sX,+ fX +*% H-:, .. 

= £ -f 2 £ 3 ^ • - * 

with A^ and given - be Eq. ; (3.14), and - " 


(B. 


^2 ^ - bn pr -Xl 


i z 


(B. 


■fj = { ‘ (2;b>np^(x« lR Xz,^ r X 2 ,p , 

+ 2 C npv ..X. lP -X^Xy- ] 

Separating terms. of the same order of magnitude .yields > 


(B. 


H- 


2 ) 


3)\ 


4 a) 


4b) 
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9X» 


A 0 X 2 =-F, 


d U 


+ A x =X - dXi _ d X i - A v 
di 6 3 9^, c>i, 

B.2 Second Order Solution 

The solution of Eq..(B.5a) is given by Eq. (3.17), 

i ^ i o ■ /' — p 

X, = UQe + C.T. 

where a is now a function of t^t^ — BY combining Eqs. 

(B. 5b) and (B.6), one obtains 
7 ) y A , l 2 oi t o - ,, p -r 

^ + A x, = q t, e + a a t - ®Ue + C. i- 

*0 l 2, 


(B. 5b) 
(B. 5c) 


(B.6) 


die 

where 


DM 


(2) 


(o) 


p'l 


ri n 


u, 



? bnpo ( U 


^ Unp^ ( Up t 


The condition to avoid secular terms is 

Qo. = o - ■ 

dl< ■ i. 

and -then -the " solution ^for _xr is given hy* — 

ur i2a>te- - _ (o) - - 

x v - Q p* e + cto-P* -+ c. V. 


(B . 7 )>. . 

{ B . &a:) zm 
(B . 8b) 

(B. 9')^, 9} 


( B . 10) ~~ 


where ” - ■ 


P 2 % [-t2a)I+ A 0 | f. 


-1 ciXZl 
Z 


(B . 11a) 


p:°-= Ao-fj.-. 


(B .lib) t 

It may - be T noted -that ' the ^second’ order 'nonlinear: -terms ..do- .-“j 
not yield ' any . secular, terms." For this reason vthe expression --- 


*Note--that; as it will .be clear., from "the : third :order-analys.is-^: - 
the term of the type ||€~ - (Solution of homogeneous -equation). 
can be safely disregarded . . . 
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used for X , Eq. (3.11), does not contain terms of order £ , 
since the perturbation of A must be of the same order of 
magnitude as the nonlinear terms which yield secular terms. 


B.3 Third Order Solution 

Combining Egs (B.4b), (B.6), and (B. 10), one obtains 

^3 { 2 (bnp^ + bn^p) X< iP Xz.c^ t C« PV - X< lP x n Xt.r } 

— f ^ / \ / *- . ~ L V x t2t Aiic 

“ | i vbnp^tbn^jVa UpC + a Up^e . A a P*<% e 

/s. In)- 2 ^ ^ / l : 

+ 2 aa a a e; ,= J + g. (a^ F e - +- 


4 

a /5^Q^-^au,e‘ 

au r e Jj 

3 tddJ'io 2 — ^ {.olio /-^ -t 

— O 4^‘ G + (X (i 7^6 4- C. 1 . 


a)4o 

4- 


where: 


^ ^ p2,c| ■•; * 5 r 0 r p^r- - U^ U r 


PV 


(B. 12 ) 


(B; 13a) 


2 thnp^.b n p)tup.p 2 !f A2 Up 


t Gnp<^<- \ Up Ur t Up U* Ur t Up Ur 


(B. 13b) 
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Finally, combining Eqs. (B.5c), (B.6), (B.9), (B.10), and 

(B.12) one obtains 

+■ 

d+o 

Lti-ie m l U)^9 


o., . i3w)^o po) _ * — 

r X . 3 +. ^ y = *f 0 6 +■ T 2 Q- ^ ^ 


-9aue^-Au^ ^ CI - 

d±* z 


(B. 13). 


which is formally identical to Eq. (3.18). [The definitions 
of and , however, are now given by Eq. (B.12) instead 

of Eq. (3.19)]. Hence, the solution of Eq. (B.13) is given in 
Section 3 ; The conclusions.- are -the one s given in Section 3r4 .i_ 
The'j solution ;is given by. iv 

where a is given in Eq, (3.24). Note that the only effect 
of the . second ...order;'no.n linear. -terms -on Eq. _(B. 174) „ is the order 
of magnitude^ dif 1 the rd errori ■: arid sttfe tdif fereht ^definition, ot cl 
f J [Eq.‘CB_-. 12) ^instead- of Eq. 113 .T9) ] ; “whichiis. used.-rin. the ir-_ 
definition -of ^ [ Eq ^ (3 . 2 3c )T- - i } - 
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APPENDIX C 

THE FUNCTION h j (t 2 ) 


C.l. Introduction 

In Section '4 , Eq.“(4.23.b) 

£ir + (3 r b, 4 / R I b, + 2 H 1 = 

<r |«| + l^l 3 + I*/ - f*/ 

was derived. - By setting,. (see Eqs. (4.24.b)) 


t ^ B h 
r r 1 


and -Using _ Eq." (4 . 25) - ; one -~ob tains- "Eq. (4;. 27 . bV - r - 

r - 2^ jtcj + S + <T 'M 1 + .Tj / W 4 
5 fc* 

C 3 f / 2fc 4 


<C.l) 


(C . 2 ) 


(C. 3) 


On the -other • hand,- • according; 'to Eq$.<4 . 36) arid (4 s 2 6). 

_ _E*. + <a 0 /-/HJS-] ] 1^1 3 

R 


■= 


R 


/*■ 




T$ M 


— I *-1 


S Hr 


dictf 
3 t 2 


cc.4 r 


while, ^according .to Eqs .:(2 ,.9),and (3-. -2 4). (see also vEq. '(2 .23)) 


® POOR 


A riTi ret 
TO 

QUA?/: k 


f = + -2L ^/a| + 

2* 


(C.5) 
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and thus 


at. 


iL JL fi-i |«f ) 


+ 


01* 


- -i*- fe-<r n Vi*\ 2 I-i ^ + ^ 




V^/ Jl*l 0fc, 


(C.6) 


£ 

since (deriving ^Jf-L from Eq. (4.32), € from Eq. (4,31) and 
from Eq. (4.35)) 


9k 


Z un - r-, JC I * 1 * -T 


■pt 4 i a l 


«i 2 € x i£ 


at. 


*r -*!*</% zKis^s; (4?) 1 ] 

-p.m**. . L w#/ 


r '^1 Jr ]i«i - 

F inallyjcombining ;Eqs ;■;(©. 3 £ G>4 OaBdsGtf^onecobtains.* 


(C.7) 


•at. 


p«* w 


j + r (&-) \i°- 

R K JK ■ R ' /* ' J 


+ 2 ?f 


B 


W 


2/«/ 


P< 




4 


<5"; * ^ + V M 



-^*1 - 




(C . 8) : 


This! is the -desired differential for IB" ^ 2 ) ;— . . This\.-_I-i 
equation ' is a n a 1-y zed < in t hi s ' Appendix v i y , 


ORIGINAL 

OP POOR 


page is 

quality 
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C;2. The functions (t) o anc * b (t^ j 
Equation (C.8) may be rewritten as 


~~ - % + 7, M 2 + % l*l 4 + % 4 1*1 




where 


^ + % -2- i*\A 


(C.9): 


1- -*■- Ik 

1. - s' . 


1, 


1 


<f* _ ?r V* 


A 


z >- 


£ - S' &+ s* 
* * 



j{r. B' 0 ' 

Z 3 * * 



C.10 



Integrating Eg. (C. 9) one obtains* 



In order to avoid secular terms, set 


l -1, h + 


- O 


(C . 12')- 


or feet B g\r(C, 10 )) . 

ORIGINAL PAGE IS 

% f - 0 OF POOR QUALITY 

9fc 4 ' ■ 


(C. 13) 


. v h) . £ £ r :, *•'*'*- 

*" a * $ ‘" w - i w ‘}= ft) 1 ^F, (fiH'*"') 
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with 


*.= [ 


<5 - S' 

i i 


Jte- 




Equation (c.13) yields 

and ^combined with^Eq. (S'. ll), z y i e Id S' " - ~ 

B = B"Vb; , V| i + 3 i w ^|«| + £f' \«\'U\*\ 


(C.14) 


(c . is); 


(C. 16) 


with .l( see rEqs 10) .arid :(4id3 71 i) 

B ( \ _ JTi * „ n 

1 - 7 z;: l 'is-n* 


OBIGINAL PAGE ® 

0F POOR QU^ Lm 


- -J. W - HLf St- - s; (£*?]} 

1 L if* * ( ?*' JJ 


4. A B c w + 1 J& < 5 .* 

PR <*R 


= .21 — l (s;'-h_ s'') 

Hr 2 S k 
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FinalXy^jCombiningj -Eqs 2kand (C;d6)-yi.elds-i«_i-L. 

I b w 1*1 + +(b;% b ( ; m 2 ) m < c - i8) 



